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Abstract. - The frequency dependence of the dielectric permitivity's maximum has been 
studied for poled and unpoled doped relaxor strontium barium niobate Sro.6iBao.39Nb206 ■ 
Cr^+ (SBN-61:Cr). In both cases the maximum found is broad and the frequency dispersion is 
strong. The present view of random fields compensation in the unpoled sample is not suitable 
for explaining this experimental result. We propose a new mechanism where the dispersion of 
quenched random electric fields, affecting the nanodomains, is minimized after poling. We test 
our proposal by numerical simulations on a random field Ising model. Results obtained are in 
agreement with the polarization's measurements presented by Granzow et al. [Phys. Rev. Lett 
92, 065701 (2004)]. 



Relaxor ferroelectrics were discovered almost fifty years ago [1] and they have recently 
found a multitude of technical applications [2]. In particular, from the point of view of 
optical devices, uniaxial strontium barium niobate (SBN) has been proved to be extremely 
useful [3,4]. Relaxor ferroelectrics present a transition with temperature from polar to non- 
polar phase which is characterized by a broad, frequency dependent maximum of the dielectric 
permitivity [5]. Another important feature of the relaxor ferroelectrics is the existence of 
locally polar nanoregions well above the Curie temperature Tc which persist until the so- 
called Burns temperature Tb [6], which is close to 350°C in SBN [7]. 

In SBN, these properties have been studied taking into account the inherent disorder 
of the relaxor [8,9] by means of a random field Ising model (RFIM) [10]. The disorder 
produces quenched local random fields which interact with local dipoles. The local random 
fields stabilize the dipoles in a small nanorcgion, leading locally to a nonzero value of the 
spontaneous polarization even above the Curie temperature. The use of a RFIM to model 
uniaxial ferroelectrics has been supported by recent findings in SBN which point toward the 
existence of internal fields [11-13]. 

Based on the RFIM, it is possible to make predictions about critical exponents in SBN. 
Theoretical calculations for the RFIM predict small values for the spontaneous polarization 
critical exponent (3 — 0.06 ± 0.07 [14]. Due to this small value, the decrease of polarization 
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when increasing temperature close to the Curie point must be very sharp compared with the 
behavior of other systems such as Ising models, where p ^ 0.325. 

Actually, there has been some controversy about the value of (3 for SBN crystals. NMR 
measurements yielded (3 = 0.14 [15] very close to the one predicted by the RFIM. However, 
linear birefringence measurements in SBN resulted in /? ~ 0.35 [16], which is actually closer 
to the value of the Ising model. This discrepancy has been recently solved by pointing out 
that NMR samples where fully poled, while linear birefringence samples where unpoled [17]. 
Granzow et al. have shown experimentally how the critical exponent /3 varies from /3 = 0.3 in 
the case of a non-poled sample to /3 = 0.126 when the sample is fully poled. This is equivalent 
to say that polarization decreases faster when increasing temperature in poled SBN than in 
non-poled SBN. 

What is the reason for this behavior? A possible answer came from Ref. [17]: Basically, 
quenched random fields are compensated by those emerging from charged fractal nanodomain 
walls in the unpoled sample. Since there is no random electric field, the behavior of the system 
resembles the one of an Ising model. However, when SBN is poled, the compensation of the 
electric fields is not so effective due to the larger size of ferroelectric domains. Since the electric 
fields do not dissapear, the behavior for the poled sample must be the one corresponding to 
a RFIM. 

Assuming that quenched random fields are responsible for the relaxor behavior in SBN, a 
sharper dielectric peak and a decrease in the frequency dispersion should be expected for the 
unpoled sample. To check this hypothesis we have performed dielectric permitivity measure- 
ments in SBN-61:Cr for poled and unpoled samples. 

The SBN-61:Cr unpoled samples were small platelets of thickness O.Smm < d < 1.8mm 
and areas between 12mm^ and ASmm^. They were doped with a small amount {Q.QlwtVo) of 
Cr^+, but it seems that the impurity type is more relevant from the optical than the ferro- 
electric point of view. The two faces perpendicular to the ferroelectric axis were covered with 
electro-conductive silver paint. The samples have been grown by Czocharalski technique at 
the General Physics Institute of the RAS, Laser Materials and Technologies Research Center, 
Moscow, Russia. The measurements of the dielectric constant were performed under IT^ of 
applied voltage, with a LCR meter HP 4284A whose frequency ranges from 2QHz to IMHz. 
The temperature control was made through a Unipan thermal controller type 680, provided 
with a sensor whose accuracy is O.IK and the sample actual temperature was obtained with 
a T thermocouple and a DMM Keithley 196. The poled sample was measured using the 
following procedure: Once the system reached a temperature higher than the critical temper- 
ature, we pre-polarized the sample applying a constant 2.8kV/cm DC field while freezing the 
sample until RT and then the dielectric constant measurements were performed for the same 
frequencies that in the unpoled case. 

The experimental measurements presented in Fig.l show the typical relaxor frequency and 
temperature dependence of dielectric constant, with a rounded peak and a non-well defined 
transition temperature. So actually, both samples behave in a very similar way from the 
relaxor behavior point of view. The value of the peak of dielectric constant we find for 
SBN-61:Cr is much smaller than the one found in Ref. [18] for SBN-61:Ce. This small value 
may be due to the quality of the sample or to the fact that Cr doping increases the SBN 
relaxor behavior. In any case, both materials present relaxor properties for poled and unpoled 
samples. 

It is not easy to explain the unpoled results when considering random fields compensation 
in the unpoled sample. The broad maximum and the frequency dispersion found for the non- 
poled sample imply that quenched random fields are not compensated. On the other hand, the 
broad maximum also appears for the poled sample, indicating that random fields are present 
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Fig. 1 - Dielectric permitivity vs. temperature for non-poled sample of SBN-61:Cr. Inset: poled 
sample 



too. However, there must be some differences between the fields distributions when comparing 
both samples, in order to explain the changes observed in the polarization behavior [17]. How 
to explain the dielectric permitivity experimental results and how to make them compatible 
with the ones observed for the behavior of spontaneous polarization versus temperature in 
Ref. [17]? 

In the following we will propose a model to answer this question. Local fields are due to 
positive charges either in deficit at unoccupied Ba^^ and Sr^^ sites in the open tungsten 
bronze structure [19], or in excess at Cr^+ ions replacing Nb^^ [20] . 

The fields are quenched due to the low mobility of the charge carriers below the Curie 
temperature [21] and they are almost uncorrelated due to the short-range character coming 
from a large dielectric dc susceptibility [9]. An electric field applied to polarize the SEN 
sample at temperatures below the transition point should interact with the charge carriers, 
introducing some degree of correlation between them [22], and decreasing the value of the 
dispersion of local fields. Once the electric field is turned off, the new "ordered" internal 
electric fields become quenched at room temperature. Of course, a change on the dispersion 
(T of the Gaussian distribution of fields in a RFIM does not change the critical behavior of 
the system (only the critical temperature) [14] so, a change on /3 is not expected. However, 
we should take into account that nanodomains found in relaxors at high temperatures are not 
formed by just a single dipole. It is necessary to consider that a quenched randon field affects 
not just to a single dipole but to a local group of dipoles forming a nanodomain. A change on 
the dispersion of local fields affecting complete nanodomains may change the behavior with 
temperature of the total polarization of the system. 

We propose the following hypothesis: (i) A poled sample presents a distribution of nan- 
odomain electric fields with an small dispersion. Due to this small dispersion, all the nan- 
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Fig. 2 - Distribution P{ha) of fields for a — 2, corresponding to the poled sample (circles) and for 
a = 3 corresponding to the un-poled sample (squares). 



odomains behave similarly when the critical temperature is approached and they decrease 
their polarization almost simultaneously. This commom behavior implies a sharp transition 
for the total polarization versus temperature (/3 decreases), (ii) A non-poled sample presents 
a distribution of nanodomain electric fields with a huge dispersion. When the critical tem- 
perature is approached the behavior of the nanodomains is very different from one to the 
other and the local polarizations do not decrease simultaneously. This behavior results on 
a smooth behavior of the total polarization versus temperature (/3 increases). However, in 
both cases, nanodomain electric fields are not compensated explaining the broad peak and 
frequency dependence found in experiments for both, poled and unpoled samples. 

To check this hypothesis we present Monte Carlo simulations of several random field Ising 
models, with a Gaussian distribution of fields centered in zero and with a width of the distri- 
bution given by a. We will consider different values of a, corresponding to poled and unpoled 
samples of SEN. To study the effect of a random field affecting not just to a single dipole but 
to a group of dipoles, we divide the system in several nanoregions a, assigning a single electric 
field to all dipoles belonging to the same nanoregion. To make calculation easier, we consider 
a cubic system formed by A'^^ x L'^ dipoles with periodic boundary conditions formed by small 
cubic-shaped nanoregions with dipoles. The Hamiltonian reads as follows: 

H = ~J ^ SiSj - Jh^^haSi (1) 

<ij> ct i£a 

and the distribution of quenched fields is given by, 

The distributions of fields considered in this work for the poled and the non-poled samples 
are shown in Fig. 2 
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Fig. 3 - (a) Polarization P{T) vs. temperature T for a = 2, corresponding to the poled sample 
(circles) and for cr = 3 corresponding to the un-poled sample (squares). Continuous lines represent 
a fitting to a power law. (b) ln{P{T)) vs. ln{Tc — T) for a = 2, corresponding to the poled sample 
(circles) and for u = 3 corresponding to the un-poled sample (squares). Continuous lines represent a 
linear fitting. 



For our simulations we set the values N = 4, L = 25, J = 1 and Jh = 0.02. In Fig.3a 

wc present results for the total norm ahzed polarization P = [1/iN^ x L^)] s i versus tem- 
perature for (7 = 2 (representing the poled sample) and cr = 3 (representing the unpoled 
sample). We have fitted the data to a power law P{T) = Pq{\ - T/Tc)'^ for T < T^. The 
values obtained are = 4.45 and /? = 0.17 for a = 2 and T,. = 4.44 and (3 = 0.26 for (7 = 3. 
These results for /3 compare with the ones obtained experimentally for SBN in ref. [17] for a 
poled sample (with initial spontaneous polarization P{20"C) ^ 10(^C/cm^)) and a non-poled 
sample respectively. Fig. 3b presents the power law fits in log scale. 

There is a difference between the model we are proposing and the former RFIM studied in 
detail by Newman and Barkema [14]. In the former RFIM each single dipole was influenced by 
a different random field and a change on a affected to the randomness of the system but did not 
affect to the critical behavior. In contrast, in the present model, the same random field affects 
to a group of N X N X N dipoles forming a nanoregion. These nanoregions may act almost 
independently from each other (as a collection of non-random Ising models) if the differences 
on the values of the fields are large (i.e. if a is large), changing the critical behavior from 
RFIM to Ising model. The value we find for /3 in the case of a poled sample is very similar to 
the one observed experimentally yet not exactly equal to the value (/3 < 0.06 ± 0.07) reported 
for the RFIM in Ref. [14], but it would become very similar in the special case of considering 
very small nanoregions {N = 1). 

It is also possible to understand this difference on the critical behavior of the polarization 
by considering the polarization of each single nanodomain versus the temperature. Results are 
presented in Fig. 4 for 50 different nanodomains. Note how, in both cases, the polarization for 
each nanodomain remains different from zero at T > Tc, as usual for a relaxor. However, the 
dispersion on the polarization curves for the poled a = 2 case (Fig. 4a) is smaller than the one 
we find for the unpoled a = 3 case (Fig. 4b). The critical temperature for each nanodomain 
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may be determined by considering the maximum of the derivative of the polarization in each 
nanoregion versus the temperature. A seasonal derivative with a period of ten points is shown 
in Fig. 4a and Fig. 4b. For a poled sample, all nanodomains present their critical temperatures 
in a narrow window around the critical temperature of the system. Since all nanodomains 
experiment the transition at almost the same temperature, the behavior of the spontaneous 
polarization is very sharp and the value of P very small. On the other hand, when the 
sample is non-poled, the critical temperatures of the nanodomains are distributed on a wider 
temperature range. This makes the total polarization not to decrease abruptly, obtaining a 
larger value of (3 than the one corresponding to the poled sample. Relaxor behavior is expected 
in both cases, since electric fields are never compensated. 

In Ref. [17] samples are polarized at room temperature, but it is also posible to polarize 

at high temperature and to perform field cooling (FC) [12]. In this special case, the relaxor 
exhibits an stable polarization at the cooling process with zero field applied. This experimental 
results are explained by Granzow et al. as follows: At high temperatures electric charge 
carriers arc very mobile and they arc able to travel through the crystal and compensate 
the local random electric fields not in accordance with the externally applied electric field. 
This process imprints an stable non-zero internal field capable of polarizing the sample at 
low temperatures. We have measured the critical exponents from the experimental results 
presented for the SBN relaxor at Ref. [12] and we have found a value of /3 = 0.15 for the 
heating process and a value /? = 0.25 for the cooling process. Our proposal is supported also 
by these results obtained for FC samples. From the point of view of our model, the field 
compensation coming from the travelling charges at the FC process implies a non-centered 
gaussian distribution of fields (a total value of the internal electric field different from zero). 
Also, the relaxor will show a short range ordering as the one already described for the poling 
process at room tepratures [22] , turning intto a decrease on the value of a and the value of 
/?(/? = 0.15). Once the system is heated to very hight temperatures, the internal electric 
field still remains [12], but the short range ordering dissapears {a turns broad again). Then 
,at the cooling process with no external field applied, the system will show the spontaneous 
polarization comming from a non-centered distribution of internal fields, with a broad value 
of (T, turning into a large value of /? (/3 = 0.25). Note how, in this case, it is imposible to 
explain this large /? value by a total compensation of the internal fileds (pure Ising model), 
since non-zero internal fields are needed to explain the existence of spontaneous polarization 
when cooling. 

In summary, the dielectric dispersion in relaxor SBN is similar for the non-poled and for 
the poled sample. For this reason, an explanation for the change on the critical behavior 
of polarization based on fields compensation in the unpoled sample is not appropriate. A 
mechanism capable of explaining both behaviors has been proposed. In this new mechanism 
a poled sample presents an smaller dispersion of quenched electric fields than a non-poled 
sample. Monte Carlo simulations show how this model leads to a different behavior of the local 
polarization in the nanodomains of the relaxor system. This difference on behavior produces 
a change on the criticality of the total polarization similar to the one found experimentally. 
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Fig. 4 - Polarization P{T) vs. T for (a) a = 2, corresponding to the poled sample (circles) and (b) 
a = 3, corresponding to the non-poled sample (squares). Grey lines represent polarization of each 
nanodomain (Pa) and dPa/dT vs. temperature for 50 different nanodomains. 
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